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Abstract 

We study the homogeneous turbulence in the presence of a constant average 
velocity gradient in an infinite fluid domain, with a novel finite-scale Lya- 
punov analysis, presented in a previous work dealing with the homogeneous 
isotropic turbulence. 

Here, the energy spectrum is studied introducing the spherical averaged 
pair correlation function, whereas the anisotropy caused by the velocity gra- 
dient is analyzed using the equation of the two points velocity distribution 
function which is determined through the Liouville theorem. As a result, we 
obtain the evolution equation of this velocity correlation function which is 
shown to be valid also when the fluid motion is referred with respect to a 
rotating reference frame. This equation tends to the classical von Karman- 
Howarth equation when the average velocity gradient vanishes. 

We show that, the steady energy spectrum, instead of following the Kol- 
mogorov law k _5//3 , varies as n~~ 2 . Accordingly, the structure function of the 
longitudinal velocity difference (Aii™) pa r"*' 1 exhibits the anomalous scaling 
( n pa n/2, and the integral scales of the correlation function are much smaller 
than those of the isotropic turbulence. 
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1. Introduction 

Although the Kolmogorov law E(k) ~ k~ 5 ^ 3 represents the main result of 
the isotropic turbulence, there are many experimental evidences and theoret- 
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ical arguments indicating that this is not the only sp ectrum observe d in th e 



fully develope d turbulence of in c ompr e ssible fluids (IBrissaud et al I ( 1 19 731 ). 



Moffat I fll978[). iGordienko et al I fcpQlh. iBaroud et al I (l2002h\ . 



For example, in IBrissaud et al 



(119731 ) and Moffat I (Il978h . it is shown 
through the dimensional analysis, that the energy spectrum in the presence 
of an average velocity gradient dU/dy can follow the law « k~ 7 / 3 in the 
inertial subrange. This is a particular result arising from the assumption that 
the energy spectrum is linear in dU/dy. More in general, assuming that the 
energy spectrum is proportional to (dU/dy) 13 with /3 > 0, the Buckingham 
theorem states that E(k) ~ K -5/3-2/3/3^ anc ^ different scaling exponent are 

po ssible. 

Gordienko et al I fl200lh studied the forced driving turbulence, where the 
forcing term can have various origins. The authors remarked that there are 
two dimensionless parameters, characterizing the forcing term, which influ- 
ence the shape of the energy spectrum and are responsible for the anomalous 
spectra. They showed that, in a certain interval of variation of one of these 
parameters, the spectrum follows the Kolmogorov law, whereas for an op- 
portune choice of it, the spectrum behaves like k~ 2 in the inertial subrange. 

These different scaling are caused by the shear rate which leads to the 
development of coherent fluid structures. These are streaky structures, due 

to the stretching of the vo rtex lines, which exhibit t he maxim um dimension 

along the stream direction ( Lee. Kim fc Moin ( 1990l )) . In Ref. ( Lee. Kim fc Moin 



(119901 ) ) , the authors remark that these streaky structures influence the span- 
wise correlation of streamwise velocity components, being a non monotonic 
function whic h becomes negative a t high values of the spanwise distances. 



Further, in lBaroud et al I (120021 ). the authors experimentally analyzed the 



statistics of the longitudinal velocity difference Au r in a closed cylindric tank 
which rotates around its symmetry axis at a given spin rate. The turbulence 
is generated by pumping the flow in the tank through two concentric rings 
of 120 holes each, placed at the bottom of the tank, where the source ring is 
the internal one. This generates an average radial flow that, combined with 
the spin rate, determines a Coriolis force whose magnitude varies with the 
distance from the rotation axis. As a consequence, a counterrotating flow 
and an average veloc ity gradient with respect to the tank frame is observed 
( Baroud et al I ( 20021 )). The authors found that Au r presents the anomalous 



scaling (Am") 



with Cn — in contrast with the Kolmogorov law 



(£ n ~ n/3), and that E(k) 



2 



The present work studies the homogeneous turbulence in an infinite fluid 
domain in the presence of an average velocity grad i ent V X U, using the finite- 



scale Lyapunov analysis, proposed by Ide Divitiis I (120101 ) . Ide Divitiis I (1201 if ) 



for studying the homogeneous isotropic turbulence. 

In the first section, we define the spherical part of the velocity correlation 
tensor Rij, and we derive the evolution equation for from the Navier- 
Stokes equation with V X U 7^ 0. 

To study the effect of V X U on the anisotropy and on R^, the evolution 
equation for the pair distribution function is derived from the Liouville the- 
orem, assuming that the statistical equilibrium corresponds to the condition 
of isotropic turbulence when the kinetic energy rate is equal to zero. From 
this equation, the steady velocity correlation tensor is expressed in function 
of the average velocity gradient and of the maximal finite-scale Lyapunov 
exponent and, in particular, the Boussinesq closure for the Reynolds stress 
is obtained. 

Finally, the equation for the spherical averaged longitudinal correlation 
function is determined, whose solutions depend on the average velocity gra- 
dient. Moreover, we show that this equation is still valid when the fluid 
motion is refereed with respect to a non-inertial rotating frame of reference. 
The steady solutions of this equation are numerically calculated for differ- 
ent Taylor-Scale Reynolds number and several results are presented. We 
found that E(k) pa k~ 2 in the inertial subrange, thus the statistical mo- 
ments (Am™) pa r^ n exhibit the anomalous scaling ( n pa n/2, whereas the 
integral scales of the longitudinal correlation function are much lesser than 
those of the isotropic turbulence. In the case of homogeneous turbulence in 
the presence of a steady shear rate, the spanwise correlation function of the 
streamwise velocity component is also calculated. 



2. Analysis 

This section analyzes the homogeneous turbulence with an uniform aver- 
age velocity gradient V X U. 

The fluid velocity, measured in the reference frame is v = U + u, where 
U = (U x ,U y ,U z ) and u = (u x ,u y ,u z ) are, average and fluctuating velocity, 
respectively. The velocity correlation tensor is defined as R^ = (uiii'j), being 
Ui and u'j the velocity components of u calculated at x and x' = x + r, where 
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the brackets denote the aver age on the statistic a l ensemble of u and u', a nd 
r is the separation distance ( iKarman-Howarth I (119381 ). iBatchelor I ( 11953I )). 

In order to determine the evolution equation of Rij, w e start from the 
Navier-Stokes equations, written for the fluctuating velocity ( IBatchelor I (119531 )). 
in the points x and x' 



dui 


duiU k 


dUiU k 


duiU k 


dt 


dx k 


dx k 


dx k 


du> 




dU' jU ' k 


du' 3 U' k 


dt 


dx' k 


dx' k 


dx' k 



1 dp 

pdxi 



1 dp' uV ,2 u > 



(1) 



where p is the fluctuating pressure and U' is 

U' = U + V X U r 



(2) 



being U and V X U assigned quantities. The repeated index indicates the 
summation with respect the same index. The evolution equation of is 
determined by multiplying first and second equation by and Ui, respec- 
tively, summing the so obtained equation s, and calculating the average on 
the statistical ensemble ( BatchelorT ( 19531 ) ) : 



dR 



dt 

being 



Ta + P u + 2vV 2 Ra - 



m 

dx k 



Rkj — 



dxk 



Rik + 



dRij 
dr k 



(U k - U' k ) (3) 



— ( UiUj (u k -U k )), Pij (T) = -{— _j 



(4) 



where d{...)/dxi = —d(...)/dri and d{...) /dx^ = d{...) jdr^. Making the trace 
of Eq. we obtain the following scalar equation 



dR 

dr k 



(5) 



where Rf k is the symmetric part of Rik, and R, defined as 



R — —Ri 



(6) 



gives the turbulent kinetic energy for r=0, H = T ti provides the mechanism 
of energy cascade, and P a ir) = a r ises fr om the fluid incompressibility 
flKarman-Howarth I fll938f ). iBatchelor I fll953h l 



Equation fl5]) incudes two additional terms with respect to the homoge- 
neous isotropic turbulence. —dUi/drkRik corresponds to the kinetic energy 
production if —dUi/dr k R ik (0) > 0, whereas dRa/ dr k {U k — U' k ) represents an 
energy transfer as this latter vanishes for r = 0. 

Assuming the condition of isotropic turbulence, dUi/dx k R ki does not give 
energy production, and the flow corresponds to a dying turbulence whose 
dissipation rate is calculated with the von Karman-Howarth equation. In 
this case the energy spectrum is a function of 

On the contrary, the combined effect of the anisotropy (i.e. (u 2 x ) ^ {u 2 ), 
( u x) 7^ ( u z)> Rki(fy 7^ 0, i 7^ k) and of V X U leads to the following condition 

g*„(0)^0 (7) 

which can correspond to an energy production, where the energy spectrum 
is not more a function of k, 2 . 



At this point, we want to determine an useful scalar equation for describ- 
ing the main properties of the velocity correlation. We start decomposing 
Rij, H and U into an even function of r = |r| (here called spherical part), 
plus the remaining term: 

Rijyxi "Tyi ^ z) -Rijy) ^Rijyxi Tyi ^ z) 

H(r x , r v , r z ) = H(r) + AH(r x , r v , r z ) (8) 
U'-U = U(r) + AU(r x ,r„r z ) 
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(9) 



being 

Rij(r) = - (Rijir, 0, 0) + R^O, r, 0) + %(0, 0, r)) 

+ ^ (%(-r, 0, 0) + i2y(0, -r, 0) + i^(0, 0, -r)) , 

fl"(r) = ~ (#(r, 0, 0) + H(0, r, 0) + H(0, 0, r)) 
o 

+ 1 (#(_ r o, 0) + #(0, -r, 0) + (0, 0, -r)) , 
o 

U(r) = ~ (U(r, 0, 0) + U(0, r, 0) + U(0, 0, r)) 
o 

+ ^ (U(-r, 0, 0) + U(0, -r, 0) + U(0, 0, -r)) 
o 

where Ai?,j(0)=A/7(0) = 0. Therefore, the Fourier transform of R identifies 
the part of the energy spectrum depending upon n 2 whose integral over the 
Fourier space gives the whole turbulent kinetic energy. Moreover, U = and 
AU = U' - U. 

The laplacian of R appearing into Eq. (jSJ) is written taking into account 
Eqs. © 

V 2 ^ = ^ + -^ + V 2 A/? (10) 
or A r or 

Substituting Eqs. (fTOjl into Eq. (jSJ), we obtain the following equation 

dR dAR \ fr l Arr n (d 2 R 2dR. 1V _ , , r) ri w , 
~dt + ~ST = 2 2 + \T)r* + r~dr ' + 2 " V A * " G ~ ^ 

where 

G = ^R ki + Go, AG = ^AJfc - A - U' h )j (12) 

where G represents the spherical part of —OR/ dr k {Uk — U' k ). As R, H and 
G are even functions of r, whereas Ai?, AH and AG are not, and since Eq. 
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(TTTj) must be satisfied at each time and for all the values of r, Eq. ( ITTj) is 
satisfied when 

I f + H^ + i^l <- d3) 

<)AR ^ + 2uV 2 AR-AG (14) 




dt 



Thanks to decomposition (jSJ), Eq. ( Tl3|) is about independent from Eq. ( [Ml) , 
thus .R is governed by a decoupled scalar equation which individually de- 
scribes some of the properties of the non-isotropic turbulence related to R(r). 
As far as Eq. (fl"4l) is concerned, AR is a complicate function of r and t which 
depends on the particular problem. 

Here, the turbulence is studied using Eq. ( 1T31) alone, whereas Rij(r) will 
be determined in function of V X U, by means of a proper statistical analysis 
of the two-points velocity correlation. 



3. Two-points distribution function 

In order to obtain the pair distribution function, consider the equations 
of motion of the various continuum fluid particles 

^ = v(x fc ,t) = v fe , fc = l,2,..JV (15) 

This is an ordinary differential system, where v(x, t) varies according to the 
Navier-Stokes equations, here written for each fluid particle in the following 
form 

^ = v fc (vi,v 2 ,...,v fc ,...,Vtf), k = 1,2, ...,N (16) 

being the acceleration of the k th particle. This form of the Navier-Stokes 
equations is obtained once the pressure is eliminated through the continuity 
equation. 

The distribution function of the velocities and of the spatial coordinates 
of all the fluid particles 

F = F(t, x 1; x 2 , ...Xtv, v 1; v 2 , ...Vat) (17) 
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is defined in the phase space r = Ti x T 2 x ... x T^, where I\ = {x^} x 
{v fc }, k = 1,2, ...N, being {x fc } = {(—00,00) x (—00,00) x (—00,00)} and 
{vfc} = {(—00, 00) x (—00, 00) x (—00, 00)} the subspaces of coordinates and 
velocities described by the k th particle. This function, which vanishes on the 
boundaries of {vfc}, {k = 1, 2, ■■ ■ N), sa tisfies the Liouville theorem associated 
to Eqs. (USD and flU flNicolis I fjl995h ) 



OF N N 

-Qt + E • v * = - E • ^ ( 18 ) 

k=l k=l 

To study the correlation between two points of space x = Xi and x' = x 2 , 
the two-points distribution function associated to x and x', is considered. 
This is the reduced distribution function calculated, by definition, integrating 
F over T 3 x T 4 x ... x T N 

Fi2) = I / •••/ F fl d * 3 M (19) 

Hence, F^ obeys to the equation arising from Eq. ( TL8l) 



dt 



+ V X F (2) • v + V X 'F (2) • v = J (20) 



where 
J — - 



r 3 JTa 



r / N N \ N 

/ E v ^ F - v ^ + E v v^-^ )U dx l dv k ( 21 ) 

•* r N \k=3 k=l J k=3 



is the rate of F^ caused by the interactions between the two particles that 
simultaneously pass through x and x', and all the other fluid particles. 

The local statistical equilibrium for the system of two fluid particles is 
expressed by the condition 

J = (22) 

To determine the expression of J, the case with null rate of kinetic energy is 
first considered. We assume that the local statistical equilibrium corresponds 
to the condition of fully developed homogeneous isotropic turbulence in an 
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infinite fluid domain whose values of momentum and kinetic energy coincide 
with those of the current condition 



F {2) du 3 du 
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F^du 3 du' 3 , 



F (2) u du 3 du 



f3 



F (2) u du 3 du' 3 , 



(23) 



F (2) u • u du 3 du 



/3 



F (2) u- u du 3 du' 3 



(2) 

That is, Fq represents the fully developed isotropic turbulence and is related 
to F^ through Eqs. (123]) . Taking into account the homogeneity, the centered 

(2) (2) 

moments of _F are constants in space, thus Fq is a function of v — U and 



U' 



F (2) (v, v'; x, x') = F ^(v - U(x), v' - U(x')) 



,(2) 



(24) 



Now, J is the rate of F^ 2 \n, u') whose variations are caused by the fluc- 
tuations (u, u') 



V dt ) ~ ~, 

\ / 11 ir 



(25) 



These var iations are determ ined using the finite scale Lyapunov analysis pre- 
sented in Ide Divitiis I ( 120101 ). where for sake of convenience, J is calculated 
for u=0, u' = u' e A< - r ^. According to the theory, |u' | << |u'|, and the varia- 
tions of F<® are calculated through Fq 2 ^ and the Frobenius-Perron equation 



(INicolis I (|1995D ) 

/oo poo pOO 
/ / F (2) (u, u' + u' )6{u' - u' e Mr)t )d 3 u' 
-oo J —oo J —oo 



(26) 



F (2) (u,u' + ti'e~ A(r) ') 



where 5 is the Dirac delta and X (r) > is the maximal Lyapunov exponent 
of finite scale fide Divitiis I f |2010l )). Being |u'e- A(r) *| « |u'| 



F (2) (u,u' + ue 



' 4- f,'p- A (* 



dF, 



(2) 



0_ . fi / e -A(r)t 



du 



(27) 



and, taking into account Eqs. flU and ((27]), J = A(r)(F (2) - F^). 
In presence of a nonzero rate of kinetic energy, J reads as 

J = A(r) (f {2) - F&) - J D (28) 

where — Jo is the rate of due to the rate of kinetic energy. Now, the 
statistical equilibrium (J = 0) differs from the isotropic turbulence, being 
J D ^0 into Eq. (12gj). 

Therefore, the evolution equation of F^ is assumed to be 

^ + V X F( 2 ) • v + V X ,F( 2 ) • v' = A(r) (f o (2) - F«) - J D (29) 

Accordingly, F^ varies depending on the boundary conditions associated to 
Eq. ( 1291) and on the initial condition. Equation ( )29|) is a partial differential 
equation, where 1/A identifies the relaxation time of the system given by a 
pair of fluid particles. 



4. Effect of the average velocity gradient 

To study the effect of V X U, consider now the homogeneous turbulence in 
a steady velocity gradient, whose pair distribution function is 

F^ (v, v'; x, x') = F (2) (v - U(x) , v' - U(x')) + 4> {2) (v, v'; x, x') (30) 

where (ffi 2 ', representing the deviation from the isotropic turbulence, satisfies, 
at each instant, Eqs. (1231) 

I j ^du 3 du' 3 = 0, 

J v J v' 

(2) u du 3 du' 3 = 0, (31) 
j J 4> {2) u- u du 3 du' 3 = 0. 

J v J v' 
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and Eq. (£2) 



A0 



(2) 



Jn- 



5F Q (2) 50( 2 ) 90( 2 ) 90(2) 



9* 



+ 



9F f 



(2) 



"Up + 



9F f 



(2) 



+ 



9x1 ' Vp 



(32) 



dx„ 



(2) 

where the spatial derivatives of _F are written in function of V X U = V X 'U' 
by means of Eq. (j21 



dF, 



(2) 



dxi 



^o (2) dUj 
dvj dxk 



dF, 



(2) 



dxi 



dFfdl^ 
dv'j dxk 



(33) 



With reference to Eq. (1321) . in an infinitive fluid domain, where V X U = 0, 
Eq. (129]) admits solutions F^ representing the homogeneous turbulence that 
decays because of Jd ^ 0. For V X U ^ 0, the last term at the RHS of Eq. 
f[3"2"j) identically satisfies Eqs. (13"T|) . thus 0^ 2 ) must also satisfy the following 
set of equations 



1 
u 
u • u 



(dF^ + + a^> + Jd , du , Av , 3 ^ (| (3 1 ) 



Of 



dX T 



dx' p " 



A sufficient condition for these five equations is 

dF m ^(2) 9 0(2) 

+ — V p + -rn- V v + J D 



dt 



dx r , 



dx' p "'' 



(35) 



Assuming that Eq. ( 135]) is true, and taking into account Eq. ( 132|) . 0^ 2 ) is a 

(2) 

linear function of the velocity gradient and F 



F® = F n (2) + 



OF^dU, dFF'dUj, 

~ Vp ~T~ n . n t-A 



,(2) 



A(r) 



dvj dx p 



dv'j dx p v 



(36) 



Of course, F^ really changes starting from an arbitrary initial condition, 
therefore Eq. (1361) represents an approximation which can be considered to 
be valid far from the initial condition. 
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At this stage of the analysis, the tensor is calculated, by definition 
Rki= I [ F {2) u k u' t d 3 u dV = R ki0 



Into Eq. ( 157)) . the integrals of d/dvj(FQ 'VpUkU^) and of d/dv'^F^ 'v p Uku'^) are 

(2) (2) 

both equal to zero as they are the integrals of F v v UkU i and F VpU^ cal- 
culated over the boundaries of {v} and {v'}, where F^ identically vanishes. 
Furthermore, taking into account that U is solenoidal, R^ is 

Rki = Rkio — {^q^T Rp i0 + ~Q^~^ kq ®*j 

where R^io is the seco n d ord er velocity correlation tensor for the isotropic 
turbulence ( Batchelor ( 19531 )) 



Rm{r) = u 2 ((/ -9)™+ 9h) (39) 

being / and g — f + 1/2 r df /dr longitudinal and lateral velocity correlation 
functions, respectively. Namely, according to the present analysis, Rki is the 
sum of the isotropic correlation tensor plus the further term due to V X U. 

For r = 0, Eq. ( 138)) provides the expression of the Reynolds stresses in 
function of V X U 

<->=" 2 M(t + g)) (40) 

where A = A(0) is the maximal Lyapunov exponent. Equation ( )40l) coin- 
cides with the Boussinesq approximation, where the link between (ukiti) and 
V X U is represented by the ratio u 2 /A which identifies the eddy viscosity In 



12 



particular, the velocity standard deviations along r x , r y and r z are 

/ 2\ 2/1 ^ dll x 
{u x ) = u 1 



A 9s 



= " 2 ( 1 " XT* ) ' (41) 

2 <9£/. 



« = « a 1 



A ftz 



3 

Therefore, = 3m 2 . According to the present analysis, V X U does not 

k=l 

modify the turbulent kinetic energy with respect to the isotropic turbulence: 
It only generates a non-isotropic condition which distributes the kinetic en- 
ergy 3/2u 2 along the three spatial directions in a different fashion, depending 
upon V X U. 

Now, Rki is calculated according to Eq. (j9]) 



and R and G are now expressed in view of Eq. (142]) 

(43) 

where S is the frame invariant scalar 

s = mf&Uk + du 1 \ 2 {Sik + Ztt) ^ = 2SA > (44) 



being = \/2(dU k /dxi + dUi/dx k ), Z ki = l/2(dll k /dxi - dUi/dx k ) and 
-G(0) > 0. 

Following this analysis, V X U influences the steady value of i?^, therefore 
it is at the origin of the non-isotropy as prescribed by Eq. ( 138|) . and of the 
kinetic energy production. Moreover, being G(r) 7^ 0, V X U is also partially 
responsible for the kinetic energy transfer and for the shape of the energy 
spectrum. 
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5. Effect of the Rotating Reference Frame 



Now, we remark that Eq. (l4"3j) can be considered still valid when 3? is a 
rotating frame of reference with a constant angular velocity O = ((f^)). In 
this case, Eq. (pQ) includes the Coriolis terms — 2f2v, — 2f2v' and the cen- 
trifugal accelerations — OOx, — f 2 0x'. These latter cause an increment of 

but do not provide any contribu- 
te 



(119671 )) 



the average pressure ( IBatchelor 

tion to Eq. fl3]) because they do not depend upon the fluid velocity. 



Coriolis accelerations modify the evolution equation of Rij, but do not alter 



the analytic al structure of Eq . 
Pair 



= (IBatchelor 



AH. 



In fact, due to incompressibility, again 
( Il953l )). Moreover, since is a symmetric tensor 
and an even function of r, whereas Qij is antisymmetric, Eq. ( ITS]) maintains 
its structure. 

As the result, the RHS of the Liouville equation (fTHj) includes the addi- 
tional terms — 2V Vfc • (F Ov^.) — V Vfe -F-Of2xfc whose integrals over T 3 x r 4 ... x 
T N are equal to zero for k > 2 as they are proportional to the integrals of 
—2 (F Ovfc) — FfiOxfc calculated over the boundary of T 3 x r 4 ... x where 
F = 0. In this case, the evolution equation of F^ includes the Coriolis terms 
and the centrifugal acceleration 



dt 



+ Vx F( 2 ) ■ v + V X ,F< 2 > ■ v' - 2V V ■ (F& ft v ) 



-2V V ■ (fw nv) - v v F( 2 ) ■ nn r = x( r ) (f (2) - f< 2 >) - j d 



(45) 



In a flow with O ^ 0, we assume that there exists a wide fluid region with 
a constant V X U 7^ 0, where the turbulence is homogeneous. This zone can 
correspond to a fluid region in proximity of the rotation axis of the frame 
dBaroud et all (l2002h ). 



Substituting the pair distribution function F^ = F (2) + in Eq. (JIBJ 
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we obtain 



+2 ^ p ( Fi2)n ^) + 2 ^7 (F {2) ^<) + ^7 (F (2) V ps n sq r q ) (46) 

dt + dt + dx p Vp + dx' p Vp 

where represents the deviation from the isotropic turbulence. 

As before, the term proportional to V X U provides null contribution into 
Eq. (123]) . so also the other ones proportional to fi, as these latter are surface 
integrals calculated over the boundary of the velocity space. Therefore Eq. 
( 134)) still holds in this case and Eq. ( 135]) is again considered to be satisfied 
for solving Eq. (1M|) . Thus, in view of Eq. gSP, 0( 2 ) is the sum of a linear 
term of V X U (as in Eq. (13"B]) ) plus another one depending on fl 

where 

^ = M^sfv + ^ ^ + ^ ^ (4g) 

Note that, <// 2 ) really varies starting from an arbitrary initial condition, there- 
fore Eq. ([4T|) - (|48p represents an approximations which can be considered to 
be valid far from the initial condition. 

The increment of the velocity correlation tensor associated to (f) n is then 
calculated 

Reiki =11 d 3 u d\> = [ [ ^u k u> d?u d\' 



2 jrp 



^ Jv Jv> dv' p 



J J (F {2) v P u k u[) - F« A {v p u k u'^ d\ d\> (49) 
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Many terms appearing at the RHS of Eq. (I49p vanish for several reasons. 
The first addend can be reduced to the sum of an integral calculated over 
the boundaries of {v} and {v'} plus the term proportional to the average 
fluctuating velocity, therefore it is identically equal to zero. The integrals 
of d/dvj(F^v p Uku' i ) and of d/dv'JF^v'ukU^) are both equal to zero be- 
cause they are the integrals of F^VpUkU^ and F^'VpUkU^ calculated over the 
boundaries of {v} and {v'}, where F^ = 0. In conclusion, Rqu is 

2 

Rflki — — ^ (^kpRpi + ^iqRkq) (50) 

Therefore 

Rq = 0, G = — - (^Rki + Rnki^j + Go (51) 
Hence, taking into account Eq. (jSTfl) 

g = - (V + (1 - (n kp s pt + n ip s pk )j + 6 (52) 

Being and Zij symmetric and antisymmetric tensors, we obtain 

dUi 

iyikpSpi + QipSpk) = (53) 



dx k 



that is, dUi/dxkRnki = and Eqs. ({13"]) and f[T3"j) are recovered. This result 
is consistent with the fact that the Coriolis acceleration is identically normal 
to the fluid velocity, thus this does not produce neither work nor energy 
transfer. 

Observe that, although R does not depend on ft, the correlation tensor 
Rij is related to O through Eq. ( )50l) . 



6. Equation for the longitudinal correlation function 

Here, the evolution equation for the longitudinal correlation function as- 
sociated to Riko, is derived. 
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Substituting Eqs. flU and flSJ into Eq. jT3), we have 



d_ 

dt 



u 

T 



or 



H „ 2 - S u 2 



3- — 

A <9r 



(54) 



+ 



Su 2 ^d\ 



3A 2 J <9r 



Grr 



First and second terms at the RHS give the mechanism of the energy cascade 
and the energy dissipation, respectively, whereas the third term, proportional 
to S, is greater than zero for r = 0, thus it provides the turbulent energy 
production. The fourth one, also proportional to S, vanishes for r = and 
does not contribute to the kinetic energy rate. 

Observe that, into Eq. ( 154"|) . H{r) is expressed without lack of generality 

as 



dK 

H{r)=3K + r— + cp 
or 



(55) 



where K, responsible for the energy cascade in isotropic turbulence, is related 
to the triple correlation funct ion k{r) = ^(x)w r (x + r))/w 3 through the 
relationship (IBatchelor I ( 119531 )) 



K{r) 



u 



d 4 
dr r 



k(r) 



(56) 



whereas ip, representing the deviation from the isotropic turbulence, is an 
assigned function of V X U and f. Ac cording to the Lyapunov theory of finite 
scale proposed by Ide Divitiis I ( 120101 ). K and A are both in terms of / 



A 



K 



u 

r 



u 



'2(1-/) 



1-/0/ 



(57) 



2 dr 

As the result, Eq. (154"]) admits the following first integral 
8u 2 f T ^ n 2 fd 2 f Adf\ 2Su 2 f 



0t 



K + 2uu z 



3 A 



+- 



, fr 3 dX 



(58) 
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being <p = ip — 2G . The first three terms appearing at the RHS of Eq. (155]) . 
identically satisfy the continuity equation ( IBatchelor I ( 119531 )) for an arbitrary 
correlation function verifying the incompressibility condition. From Eq. 
and assuming that the integral of r 2 (p converges, the last term of Eq. 
behaves like r~ 3 for r — > oo. Therefore, the continuity equation is satisfied 
when 



2 2 fr 3 dX 2 
-Su 2 



A 2 dr 



+ r (p 1 dr — 



(59) 



at each time and for arbitrary / satisfying the continuity equation. Being (p 
depending upon / and V X U, the integrand of Eq. (1591 identically vanishes 
and is in terms of / and S 



l S u 2f —— 
3 A 2 dr 



(60) 



As the consequence, Eq. ( 1581) becomes 



du 2 f 



K + 2uu 2 



9V + 4a/ 
dr 2 r dr 



2 Su 2 f 
3 A 



(61) 



whose boundary conditions are ( iKarman-Howarth I (1l938l ). IBatchelor I (1l953l )) 

/(0) = 1, lim/(r)=0 (62) 



For V X U —7-0, Eq. (1611) tends to the classical von Karman-Howarth equation. 

Putting r = into Eq. ( 16T1) . we obtain the equation of evolution of the 
turbulent kinetic energy 



du 2 
~dT 



s 



u 



(63) 



which states that 



S ^ 15 , du 2 
— — — then — — —0, respectively 
A 2 <R T dt < 1 y 



(64) 



being A^ = a/— 1//o and Rt = u\t/v the Taylor scale and the Taylor scale 
Reynolds number, and A = u/Xt- 
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The Fourier Transforms of / and K are ( Batchelor ( 19531 )) 



■ E{k) - 


i r°° 


- u 2 f{r) - 




Jo 


K(r) 



2 2 

K r 



sin kt 



ft?" 



cos kt dr 



(65) 



where E(k) and T(k) are the parts of the energy spectrum and of the so 
called transfer function which depend on n 2 . 

This analysis states that the turbulence is described by the isotropic cor- 
relation tensor R^o (i.e. /(r)), and by a non-isotropic term which depends 
upon V X U. The determination of R^. has been reduced into the following 
two steps: 1) Calculation of / and Rn-a through Eq. ( l6T|) in function of S. 
2) Calculation of Rn~ in terms of V X U and fl with Eqs. ( 1381) and ( 1501) . 



7. Steady Solutions 

According to Eq. 



(66) 



the steady correlation functions are obtained for 

S 15 
A2 ~ R^ 

with df /dt = 0, thus the dimensionless equation of / is 

'd 2 f , 4df 



l-f df 2 
2 rff i? T 



df 2 ^ f df 



+ 



10 



/ r 



o 



(67) 



where f = t/Xt and the boundary conditions are expressed by Eqs. (|62|) . As 
the solutions f E C 2 [0, oo) tend to zero for r — >• oo, with d 2 f(0)/df 2 = — 1, 
the boundary conditions (|62|) can be substituted with the condition 



/(0) 



g(g) 
df 







(68) 



Therefore, the boundary problem represented by Eqs. ( 1671) . ( 1621) . is replaced 
by the following initial condition problem 



df_ _ p 

df 



dF 

df 



5/r 
V / 2(W) 




(69) 
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whose initial condition is 

/(0) = 1, F(0) = 



(70) 



8. Results and Discussion 



It is wo rth to remark tha t Eq. (!67j) is quite similar to the equation 
obtained by lde Divitiis I (120111 ) in the case of self-similar isotropic turbulence. 
There, the equation was 



1 ^ + JL 

df R T 



df 2 f df 



Rt 







(71) 



where again, d 2 f(0)/df 2 = —1 and the boundary conditions are expressed by 
Eq. (j68p . Equation fl671) differs from Eq. flTTj) by the presence of the last term, 
that, into Eg. (I7T|) . arises from du 2 /dt under the self-similarity hypothesis. 
de Divitiis I (J201l|) shows that /-Ik r 2 / 3 , where 2/i? T (d 2 f/df 2 + A/fdf/df) 
is negligible with respect to the other terms (see Eq. ( 171 1) 1 . This deter- 
mines that E(k) 



-5/3 



((Ait r ) n ) 
ofEq. (E 



in the inertial subrange and that the scaling law 
is satisfied for ( n ~ n/3. Here, the presence of third term 
modifies the correlation function, resulting / - 1 « r and this 
implies that E(n) ~ k" 2 in the inertial subrange and that ( n ~ n/2. 

To study in detail the energy spectrum, several numerical solutions of 
Eqs. ( )69|) are calculated for different values of Rt by means of the fourth- 
order Runge-Kutta method. The step of integration is chosen on the basis of 
the behavior of Eqs. (|69|) at the small scales when r — > oo. There, d 2 f/df 2 
is about equal to K 



df 2 



Rt df 
2^2 df 



(72) 



This equation suggests t hat Af = v2/ R can be an adequate step of integra- 
tion (jHildebrand I (119871 )) for the accuracy of the numerical solutions of Eq. 

W- - 

As in Ide Divitiis I (j201l| ). the cases here studied correspond to Rt = 100, 
200, 300, 400, 500 and 600. 

Figures [1] and |2] show double and triple longitudinal correlation functions. 
Due to the mechanism of energy cascade expressed by K ( see Eq. (J57|) 1. 
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the tail of / rises with in agreement with Eq. (I58p and this determines 
that, for an assigned value of At, all the integral scales of f increase w i th Rj . 



Comparing these results with the corresponding data of Ide Divitiis I ( 1201 lh . 



we found that, the integral scales here calculated are about one- sixth those 
obtained in the case of homogeneous isotropic turbulence. This means that 
V X U contrasts the mechanism of the energy cascade, making this mechanism 
more limited in scales. 

According to Eq. (|57|) . k decays more slowly than / and its characteristic 
scales increase with R? as prescribed by Eq. ( |69|) . The maximum of \k\ is 
about 0.06 and t his is in very goo d agreement with the numerous data of 



the literature (see iBatchelor I ( 119531 ) and Refs. therein ). Again, th e scale s of 



variation of k are almost one-sixth those calculated in Ide Divitiis I ( 1201 if ). 

Figures |3] and H] show E(k) and T(k) calculated with Eq. ( )65|) . Because 
of the fluid incompressibility (see Eq. (159]) ). E(k) w k 4 near the origin, 
whereas E(k) « k~ 2 in the inertial subrange in contrast with Kolmogorov 
law k -5//3 . This disagreement, caused by V X U, makes the mechanism of 
energy cascade weaker with respect to the isotropic turbulence, in agreement 
with the previous observation, determining a higher absolute slope of E(k) 
in the inertial range. For what concerns T(n), as K does not modify the 
fluid kinetic energy, J^°T(K)dK = (see Fig. @|. 

Note that these cases can represent homogeneous turbulent flows with 
uniform steady shear rate 



V X U 



dU 
dy 



1 





(73) 



where x and z are streamwise and spanwise coordinates, respectively. In 
this situation we expect that V X U leads to the development of coherent 
structure s, similar to streak s , cau sed by the stretching of the vortex lines 



along x (ILee. Kim fc Mo in I (Il990h l This influences the spanwise correla- 
tion function of the streamwise velocity Rn(r z ) which is here calculated 
through Eq. ( 13 8 p . once / is known. The results, shown in Fig. give 
R n (r z ) in terms of the dimensionless coordinate r* = r z /(v / 'y/S) 1 ^ 2 and 
show that R\i{r z ) intersects the horizontal axis and remains negative for 
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r z — > oo. These negative values imply a wide distribution of spacings be- 
tween the different streak y structures, whose mean value depends on R T 



( ILee. Kim &: Moin I (I1990I )). From the figure, this average distance results to 
be about 70 -f- 150 in units of viscous sca le (is/y/S) 1 / 2 , comparable with the 
results of Ref. f lLee. Kim fc Moin I fll99p[ )l 



Figures |6]and [7] illustrate skewness and flatness of Au r for the same values 
of Rt- The skewness H 3 is first calculated with Eq. flTHj) and the flatness 
has been thereafter determined usin g Eq. (T75D. Alth ough -£^(0) does not 
depend upon the Reynolds number ( de Divitiis ( 2010[ )) (see Appendix), for 



r > At, Hs{f) rises with Rt and goes to zero for r — > 00. The constancy 
of -^3(0) and the quadratic terms of Eq. (ITT)) cause that the intermittency 
of Au r increases with Rt, whereas the spatial variations of H 3 and are 
the result of fir), k{r) and of Eq. ( 17T|) . thus also the scales of variation of 
skewness and flatness are significantly smaller in comparison with those of 
the isotropic turbulence. The quantity H4 — 3 is significantly greater than 
zero for r = and tends to zero as r — > 00 more rapidly than H 3 . 

The sp atial structure of A u r is studied with the analysis reported in the 
Appendix fide Divitiis I t01$ )). using the previous results. The statistical 



moments of Au r (r) are expressed in function of r through different scaling 
exponents ( n 

((Au r ) n ) = A n r^ (74) 

In order to calculate ( n , (Ait™) are first calculated in function of f, through 
Eqs. ( JTBj) (see for instance Fig. |8]). The scaling exponents are identified 
through a best fitting procedure, in the intervals (fx, f 2 ), where the endpoints 
f ! and f 2 have to be determined. The calculation of ( n and A n is carried out 
through a minimum square method which, for each statistical moment, is 
applied to the following optimization problem 

J n (Cn,A n )= [\((Au r ) n }-A n r^) 2 dr = mm } n = l,2,... (75) 



where ((Au r ) n ) are calculated with Eqs. flT8|) . fi is assumed to be equal to 
0.1, whereas f 2 is taken in such a way that C2 = 1 for all the values of Rt- The 
so obtained scaling exponents, shown in Table [TJ are co mpared in Fig. M (soli d 



symbols) with those of the Kolmogorov theories K41 ( Kolmogorov ( 194l[ )) 
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Table 1: Scaling exponents of the longitudinal velocity difference for several Taylor-Scale 
Reynolds number. 

R T 100 200 300 400 500 600 



Ci 


0.49 


0.50 


0.50 


0.50 


0.50 


0.50 




1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


Cs 


1.42 


1.41 


1.41 


1.41 


1.41 


1.41 


Ci 


1.84 


1.83 


1.82 


1.82 


1.83 


1.83 


c 5 


2.22 


2.19 


2.19 


2.19 


2.20 


2.20 


Ce 


2.59 


2.56 


2.56 


2.56 


2.57 


2.58 


Cr 


2.94 


2.89 


2.90 


2.90 


2.91 


2.93 


Cs 


3.28 


3.23 


3.23 


3.24 


3.26 


3.29 


Co 


3.60 


3.54 


3.55 


3.56 


3.59 


3.63 


ClO 


3.92 


3.85 


3.87 


3.88 


3.93 


3.98 


Cn 


4.22 


4.15 


4.17 


4.19 


4.26 


4.32 


Cl2 


4.53 


4.44 


4.47 


4.51 


4.58 


4.66 


Cl3 


4.81 


4.72 


4.75 


4.81 


4.91 


5.00 


Cm 


5.10 


5.00 


5.04 


5.11 


5.24 


5.35 


Cl5 


5.39 


5.28 


5.33 


5.42 


5.56 


5.69 



(dashed line) and K62 ( Kolmogorov ( 19621 )') (dotted line), and with the 
exponents calculated by She-Leveque ( 19941 ) (continuous curve). We found 
that, near the origin ( n ~ n/2 instead of the expected law ( n ~ n/3, whereas 
for n > 4, ( n behaves like a m uftiscaling exponent. This is in agreement 



with the experimental results of iBaroud et al I (120021 ). where the turbulence 



is produced in a cylindrical rotating tank. The multiscaling behavior is the 
consequence of the combined effect of the quadratic terms of Eq. ( 177|) and 
of the mechanism of the energy cascade expressed by K through Eqs. (I5T 



Next, to analyse the statistics of du r /df, the distribution functions of 
du r /df are calculated using Eqs. (!83|) and (ITT)) . These are obtained with 
sequences of the variables £, rj and C generated by gaussian random numbers 
generators. The distribution functions are then calculated through the sta- 
tistical elaboration of these data and Eq. (jTTJ). The results are shown in Fig. 
[TUa and [TUb in terms of the dimensionless abscissa 

du r jdf 



{(dur/dfyy/ 2 
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where the dashed curve represents the gaussian PDF. These distribution 
functions are normalized, in order that their standard deviations are equal 
to the unity. In particular, Fig. [TOb shows the enlarged region of Fig. ITUa . 
where 5 < s < 8. The tails of the PDFs change with Rt in such a way that 
the intermittency rises with Rt according to Eq. f!77|) . 



9. Conclusions 



The equation of the steady spherical correlation function is obtained in 
case of homogeneous turbulence in the presence of an uniform average ve- 
locity gradient, and an expression of the velocity correlation tensor in terms 
of this gradient is derived. When r = 0, this tensor recovers the Boussinesq 
closure of the Reynolds stress. The solutions of this equation, calculated with 
the closure based on a specific Lyapunov theory previously proposed, allows 
to determine the statistics of Au r . In particular: 

• The energy spectrum, satisfying the continuity equation, follows the 
law k~ 2 in the inertial subrange whose size increases with the Reynolds 
number. This contrasts with the Kolmogorov law and represents a 
more tenuous mechanism of the energy cascade. Accordingly, the scal- 
ing exponents of the moments of velocity difference vary according to 
( n m n/2 instead of ( n ~ n/3, and for n > 4, these exponents exhibit 
multiscaling behavior. 

• The effect of the average velocity gradient, going against the mecha- 
nism of energy cascade, makes the integral scales of / and the scales 
of variations of k, significantly smaller than those calculated for the 
isotropic turbulence. 

• In case of uniform shear rate, the spanwise correlation function of the 
streamwise velocity component, exhibits the typical shape caused by 
the streaky coherent structures due to the vortex stretching. 

T hese results, which repres ent a furthe r appli cation of the analysis presented 
de Divitiis I (I2OIQI) and [de Divitiis I jgOllf ). are comparable with the di- 
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rect simulations of iLee. Kim Sz M oin I fll990h and are in agreement with the 



theoretical a rguments of iGordienko et al I ( 120011 ) and with the experiments 



presented in iBaroud et al I (120021 ) . 



24 



10. Appendix: Statistics of velocity difference 



In this appendix, we recall the main results of Ide Divitiis I (120101 ) . For 
sake of simplicity, we only consider the statistics of the longitudinal velocity 
difference Au r associated to Riko (or /) and to Eq. ( joTT) . This approximation 
allows to calculate all the dimensionless statistical moments of Au r once 
known the s kewn ess H 3 (r). This latter is calculated in terms of / and k 
dBatchelor I (Il953f n 



H 3 (r) 



6k(r) 



((Au r ) 3 } 

((Aw r ) 2 ) 3 ^ 2 (2(1 - f(r))f 2 



(76) 



hence , H 3 (0) = — 3/7 does not depend upon the Reynolds number (Ide Divitiis 



(120101 )). and the higher order moments are consequentely determined, taking 
into account that Au r is analytically expressed as 



Au r . 



e + ^(x(r/ 2 -i)-(c 2 -i)) 

1 + 2^ 2 (1 + x 2 ) 



(77) 



Equation f J77|) . arising from statistical considerations about the Fourier-transformed 
Navier-Stokes equations, expresses the internal structure of the fully devel- 
oped turbulence, where £, r] and ( are independent centered random variables, 
each distributed following the gaussian distribution functions p{rj) and 
p(() whose standard deviation is equal to the unity. The moments of Au r 
are easily calculated from Eq. (1771) 



E 

fc=0 



((A^) w ) 1 

{(Au r f) n/2 (l + 2^(i + x 2))„/2 



(78) 



where 



((x(^ 2 -i)-(c 2 -i)) fc ) = E( ■ ) 



(-xr<(c a -i)w -i)*-*> 



(79) 



(-l)V (i ~°> 
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In particular, H 3 , related to the energy cascade, is 



8^ 3 (y 3 ~ 1) 
(1 + 2^ 2 (1 + X 2 )) 



2^ 3 / 2 



(80) 



and ip is a function of r and of the Reynolds number Ide Divitiis I ( 120101 ) 



^(r, R 7 



Ri 



15V15 



(81) 



being ip(r) determined through Eq. ( 1801) as soon as Hs(r) is known. The 
parameter x is a l so a function of Rt and is implicitly calculated putting 
r = into Eq. (IBOI 



8^o 3 (l-X 3 ) 



(1 + 2^ 2 (1 + X 2 )) 



3/2 



3 
7 



(82) 



where ^ = V(^t, 0) and = 1.075 Ide Divitiis I (120101 ). 

The PDF of Au r can be formally expressed through t he gauss i an P DFs 
P(0' P(^) an d P(0? using the Frobenius- Perron equation (iNicolis I ( 119951 )) 



F(A«' P )= / / fp(Z)p(r))p(0 6(Aur-Au'r) d^d V d( (83) 

</£ </r/ 
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Figure 5: Spanwise correlation function of the streamwise component of velocity for several 
Taylor-Scale Reynolds numbers. 
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Figure 6: Skewness of Au r at different Taylor-Scale Reynolds numbers. 
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Figure 9: Scaling exponents of Au r fo r several Rt- Solid symbols are for the present data. 
Dashed li ne is for iKolmogorov ( 1941 ). Dotted line is for Kolmogorov (1962). Continuous 
line is for IShe-Leveaue I (|l994f l 
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